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Abstract 

Finite-size scaling, finite-size corrections, and boundary effects for critical 
systems have attracted much attention in recent years. Here we derive ex- 
act finite-size corrections for the free energy F and the specific heat C of 
the critical ferromagnetic Ising model on the M x 2N square lattice with 
Brascamp-Kunz (BK) boundary conditions [J. Math. Phys. 15, 66 (1974)] 
and compared such results with those under toroidal boundary conditions. 
When the ratio £/2 = (M + Vj/lN is smaller than 1 the behaviors of finite-size 
corrections for C are quite different for BK and toroidal boundary conditions; 
when ln(£/2) is larger than 3, finite-size corrections for C in two boundary 
conditions approach the same values. In the limit N — > oo we obtain the 
expansion of the free energy for infinitely long strip with BK boundary con- 
ditions. Our results are consistent with the conformal field theory prediction 
for the mixed boundary conditions by Cardy [Nucl. Phys. B 275, 200 (1986)] 
although the definitions of boundary conditions in two cases are different in 
one side of the long strip. 

PACS numbers: 05.50+q, 75.10-b 
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I. INTRODUCTION 



In the study of phase transitions and critical phenomena, it is extremely important to 
understand finite-size corrections to thermodynamical quantities. In experiments and in 
numerical studies of model systems, it is essential to take into account finite size effects in 
order to extract correct infinite-volume predictions from the data. Finite-size scaling [l|-|3| 
concerns the critical behavior of systems in which one or more directions are finite, even 
though microscopically large and is valuable in the analysis of experimental and numerical 
data in many situations, for example, for films of finite thickness. As soon as one has a 
finite system one must consider the question of boundary conditions on the outer surfaces 
or "walls" of the system. The systems under various boundary conditions have the same 
per-site free energy, internal energy, specific heat, etc, in the bulk limit, whereas the finite 
size corrections are different. To understand the effects of boundary conditions on finite-size 
scaling and finite-size corrections, it is valuable to study model systems, such as percola- 
tion model [||] and the Ising model 0-0. Therefore, in recent decades there have many 
investigations on finite-size scaling, finite-size corrections, and boundary effects for critical 
model systems JTIJ-2I]. Of particular importance in such studies are exact results where the 



analysis can be carried out without numerical errors. 

The Ising model has exact solutions on finite lattices with many kinds of boundary 
conditions, including cylindrical ||, toroidal and Mobius strip and Klein bottle fL6f . 



This class also includes the special boundary conditions introduced by Brascamp and Kunz 
0. The calculation of the exact partition function of the two-dimensional Ising model in 
the zero field wrapped on the cylinder was performed by Onsager in 1944 ||. Exploiting 
the exactly known partition function of the two dimensional Ising model on finite M x N 
square lattice with toroidal boundary conditions 0, Ferdinand and Fisher computed the 
finite-size corrections to the free energy, the internal energy, and the specific heat up to 
order N~ l . Recently, there has been much effort in understanding the behavior of finite-size 
corrections of the free energy, internal energy and specific heat. Izmailian and Hu 1TH] and 



Salas [17] extended the results of |7fl for the free energy and the internal energy up to order 



iV -5 and for the specific heat up to order iV~ 3 . Lu and Wu [|TBJ obtained expressions for the 
partition function of the Ising model on an quadratic lattice embedded on a Mobius strip 
and a Klein bottle. They find finite-size corrections for free energy to order N^ 1 . Brascamp 
and Kunz || calculated the partition function of the Ising model on the M x 2N square 
lattice for special boundary conditions shown in Fig. 1. Recently Janke and Kenna [|H| has 
calculated the finite-size corrections of the specific heat for this boundary condition up to 
M~ 3 order. Very recently, Ivashkevich, Izmailian and Hu |2(J provided a systematic method 
to compute finite-size corrections to the partition function and their derivatives of the Ising 
model on torus. Their approach is based on an intimate relation between the terms of the 
asymptotic expansion and the so-called Kronecker's double series [J2l| which are directly 
related to elliptic theta functions. Expressing the final result in terms of theta functions 
avoids messy sums (as in some earlier works) and greatly simplifies the task of verifying the 
behavior of the different terms in the asymptotic expansion under duality transformation 
M <-> iV '. Using this approach, Salas [[H[ computed the finite-size corrections to the free 
energy, internal energy and specific heat of the critical Ising model on a triangular and 
honeycomb lattices wrapped on a torus. 
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Using the exact partition of Ref. || and the method of Ref. pOfl , in the present paper 
we derive exact finite-size corrections for the free energy F and and the specific heat C of 
the critical ferromagnetic Ising model on the A4 x 27V square lattice with Brascamp-Kunz 
(BK) boundary conditions || and compared such results with those under toroidal boundary 
conditions |7|JT5|. We find that when the ratio £/2 = (Ai + 1)/2M is smaller than 1 the 
behaviors of finite-size corrections for C are quite different for BK and toroidal boundary 
conditions; when ln(£/2) is larger than 3, finite-size corrections for C in two boundary 
conditions approach the same values. In the limit M —>■ oo we obtain the expansion of the 
free energy for infinitely long strip with BK boundary conditions. Our results are consistent 



with the conformal field theory prediction for the mixed boundary conditions by Cardy fTl 
although the definitions of boundary conditions in two cases are different in one side of the 
long strip. 

This paper is organized as follows. In Sec. II we show how to lead the partition function 
of Ising model under Brascamp-Kunz (BK) boundary conditions to the form of partition 
function with twisted boundary conditions. In Sec. Ill asymptotic expansions of the free 
energy is presented. In Sec. IV expansions of the specific heat is presented. Our results are 
summarized and discussed in Sec. V. 



II. ISING MODEL UNDER BRASCAMP-KUNZ 
BOUNDARY CONDITIONS 

For the Ising model on a lattice G of N sites, the i-th site of the lattice for 1 < i < N is 
assigned a classical spin variable Sj, which has values ±1. The spins interact according to 
the Hamiltonian 

(3H = -JJ2 SiSj (1) 

<ij> 

where J is exchange energy, the sum runs over the nearest neighbor pairs of spins, and 
(3 = 1/ksT is the inverse temperature. The partition function of the Ising model is given 
by the sum over all spin configurations on the lattice 

z Ismg {j) = Y.^ H(s) - (2) 

s 

As is mentioned in introduction there are a few boundary conditions for which the Ising 
model has been solved exactly. Among them is the special boundary conditions studied by 
Brascamp and Kunz (BK) |§. They considered a lattice with 2M sites in the x direction and 
Ai sites in the y direction. The boundary conditions are periodic in the x direction; in the 
y direction, the spins are up (+1) along the upper border of the resulting cylinder and have 
the alternative values along the lower border of the resulting cylinder as is shown in Fig. 1. 
Recently Janke and Kenna |19[ has analyzed the Ising model in two dimensions with these 



boundary conditions. They have derived exact expressions for the finite-size scaling of the 
specific heat up to the Ai~ 3 order. In this paper we obtain all exact finite-size corrections 
for the free energy and the specific heat. Moreover in our case the terms of asymptotic 



expansion are analytical functions. They are related to Kronecker's double series |22[ which 
in turn can be expressed by elliptic ^-functions [E(J. 
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For the BK boundary conditions, the Ising partition function given in Ref. can be 
rewritten as 

M M 

wH^ef^nn^j) (3) 



where \i = l/21nsinh2J and 



F(i,j) = A 



O • i 2 , • 2 ( 7T(i- 1/2)^ 2 / 7Tj 

2 smh /i + sin — — = + sin 1 



2Af J \2{M + 1) 



(4) 



Now we try to express the partition function Z M ^M given by Eq. (||) to the form of partition 
function with twisted boundary conditions Z a} p(fj) 



N-l M-l 



= n n 4 [ sin2 Br 1 ) + sin2 (^) + 2 sinh ^ ] (5) 

i=0 j=0 



for which a general theory about its asymptotic expansion has been given in Ref. ||20|| . For 
this purpose we can express the double products 

n 2 f _1 n^ +1 F(i + i, j) through nf =1 n£i F(i, 3 ) 



as 



2JV-1 2M+1 ( M M \ 4 2AT-1 

n n^+u)= nn« n^+^i^+i^+i)- (6) 

i=0 j=0 \i=lj=l J i=Q 

Here we use the properties of function F(i,j) 

F(W+l-k,j) = F(k,j) and F(i, 2M + 2 - k) = F(i, k). (7) 

This transformation leads the rectangular lattice Ai x 2J\f under consideration to the lattice 
2(Ai + 1) x 2j\f. In what follows we will use for convenience the definition of the aspect 
ratio as £ = =^M- instead of conventional one (£ = ^=). 

The left-hand side of Eq. (H) is nothing but the partition function with twisted boundary 
conditions Zf /2fi (jJL) given by Eq. fl|) with N = 2J\T and M = 2(.M + 1). With the help of 
the identity J2§ 

J] 4[sinh 2 w + sin 2 (^± 1 ^)] = 4cosh 2 (2AM 

i=0 

the second product in the right-hand side of Eq. (Q) can be transformed into the form 

2M-1 

IJ F(i + 1, 0)F(z + 1,M + 1) = [4 cosh (27V^(0)) cosh (2A/"^(tt/2))] 2 (8) 

i=0 

where 



w M (A;) = arcsinhysin fc + 2sh (9) 

is a lattice dispersion relation. 

By the aid of Eqs. (f|) - @ and (§) the partition function Z M ^M can be expressed as 
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2 (^2e^J 
Zm ^ = 4cosh[2AT^(0)]cosh[2Ar^(7r/2)] Z W^' ^ 
For our further purposes we transform the partition function ^1/2,0 m to the simpler form 

JV-l 

Z 1/2 ,M = II 2sinh [%(^)] (11) 



n=0 

where iV = 2JV and M = 2(M + 1). 

III. ASYMPTOTIC EXPANSION OF THE FREE ENERGY 

In the previous section it was shown that the partition function of the M. x 2J\f Ising 
model with BK boundary conditions can be expressed in terms of the partition function 
with twisted boundary conditions Zu 2l o, which has been well studied in Ref. [p0| . Farther 



we will use it and for simplicity we will remind some necessary parts from there. For 
reader's convenience, all the technical details of our calculations and definitions of the special 
functions are summarized in the appendices at the end of this paper. Considering the 
logarithm of the partition function with twisted boundary conditions, Eq. ([TT|) , we note, 
that it can be transformed as 

InZ 1AO (0) = ME^ffl + X>(l - e- 2 ^(^)) . (12) 

n=0 n=0 ^ ' 

The second sum here vanishes in the limit M — > 00 when our lattice turns into infinitely long 
cylinder of circumference N. Therefore, the first sum gives the logarithm of the partition 
function with twisted angle 1/2 on that cylinder. Its asymptotic expansion can be found 



with the help of the Euler-Maclaurin summation formula [24 



=li;^ - <Bf - *± (*)'^ g|. (13) 



Here S = NM = 4J\f(Ai + 1), B*/ 2 are the Bernoulli polynomials B" at a — 1/2 which are 
related to the Bernoulli numbers B p = B° p as B\l 2 = (2 1_p - 1)B P and f = § = We 
have also used the symmetry property, uo(k) = uo(tt — k), of the lattice dispersion relation 
given by Eq. (0) and its Taylor expansion 



00 \ 

A2p 



where A = 1, A2 = —2/3, A4 = 4, etc. 

We may transform the second term in Eq. fli~2] ) as 

Nl < 6. 



Eln(l-e- A M^)) = ln^ + ^B 

n=0 ^ ' V 



1/2 
2 

77 



00 A 2p Re K^°(zAO-B^ 2 
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where rj = (6 I 2 6 I 36 I 4/2) 1//3 is the Dedekind-r/ function; 2) #4 are elliptic ^-functions and 
K 2 p+ 2 (iA£) are Kronecker's double series |]20l , |22| (see also Appendix |A]). Taking into account 
the relation between moments and cumulants (Appendix |B[), the differential operators A 2p 
that have appeared here can be expressed via coefficients A 2p of the expansion of the lattice 
dispersion relation as 



A 2 
A4 

A fi 



A 2 , 

d 

A 4 + 3A^- : 



d_ 



d 2 



\q + 15A4A 2 — + 15A 2 



p h^w) '"Wj h\...k r \ d\ k - 



Here summation is over all positive numbers {ki . . . k r } and different positive numbers 
{pi, . . . , p r } such that piki + . . . + p r k r = p and k = k± + . . . + k r — 1. 



Substituting Eqs. fllTf ) and (15) into Eq. (|T2|) we finally obtain exact asymptotic expan- 
sion of the logarithm of the partition function with twisted boundary conditions in terms of 
the Kronecker's double series 



In Zt/2 o(0) — — [ uJ (x) dx + In — 

7T JO 11 



_9 t^ltiY ^ ReI Q^) 

S J (2p)\ 2p + 2 



(16) 



P =i 



where JqUJ (x) dx = 2G and G = 0.915966 is the Catalan's constant. 

After reaching this point, one can easily write down the exact asymptotic expansion of 
the free energy, F = — In Z M 2 N, a t the critical point. Plugging Eq. ([16]) back in Eq. (|T0"1) 
we have finally obtain 



-2MA/"(^ln2 + — ) +2M 

>2 7T / 

.1/2,0, 



f 7T \ 



2 p A 2p ReK 2p+2l 



ll n( l + V2)-^ 

.2 7T 



-in — 

2 2rj 



2p 



(17) 



1 o 

Note that the Kroneker's functions Kp (iA£) can be expressed in terms of the elliptic 
^-function only. Thus, Eq. ( |17|) can be rewritten in the following form 



F = 2MMf bu i k + 2JV/1 + / + E 



P tl(2yV) 2p ' 



(18) 



where 
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fbulk 

h 

fo 

h 



iln2-^ 

2 7T 



iln(l + V2)-- 

Z 7T 
1, 04 

— m — , 

2 2r?' 

vr 3 e n 



0.929695..., 

-0.142435..., 



360 V8 
48384 



/)8 MM 
#4 - U 2 U 3 J > 



5 

16 



+ {0\ + 6$) 



31 
16 



#4 + ^2^3 



l + 4£ 



4 

02 



/6 



87091200 



/ fll2fl4 
7O7r 2 £ 2 4 4 0^ + 2 4 



810S0 



2^4 



2 8 

295040 12 6350 16 \ 



4 
2 U 4 



+ 63Ott£0 4 4 e 



'2 v i 

16 



8 64 

^M/il2 '2'/il6 

— Or," A (7/1 

4 2 4 128 4 



64 ; 

' ft' s 

1+4 t 



.04 



711 + 5040f 7p + 8400£ 
02 



ail 

2^2_ 
2 



560£ 



2^2 



(19) 
(20) 
(21) 

(22) 
(23) 



(24) 



h 



^ 12 M 



C7 4 C7 3 



, 128O0 24 + 2O2240 2O 4 + 83 6 6 40 2 16 0f 
33634123776 V 2 24 24 

+ 21O49 60 2 12 12 + 3611150^0 16 + 32391O0 4 0f + 1O731O0 24 



7T i 



ejei (i + 4j|) 

1868562432 



(128O0 20 + 88320 2 16 4 + 19O560 2 L2 4 5 



+ 335680 2 *0 12 + 386550 4 0f + 153300, 



2^4 



7r lo e 2 04 4 3 4 



1177194332160 



f 0/ n/2 011 

3789 + 27720£- 2 - + 48720^-^ + 2240£ 2 2 



128O0 2 16 
vr 9 e (0 2 4 4 



02 02 02 

8/18 , r-i^cMfl io i oc:c:p;fll6^ 



31360 2 12 4 + 32 1 60*0* + 51760*0 4 12 + 25550^ 



3 4 ) 



235438866432 

,0o 



1479 + 15156£ 



$ 




0/2 0/3 



+ 2520£ 2 - 2 - + 7980^/ + 140^ 
2 



\Q_2_ 0J 

02 02 



+ 47880f 

02 



+ 47880^ 



05 



(25) 



The free energy per unit length of an infinitely long strip of width L at critically has the 



finite-size scaling form |1C 
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F = fL + r + j + ..., (26) 
where / is the bulk free energy per unit area, -/* is the surface energy, L _1 is a scaling field, 



A = 


7T 




~12 


A = 


7T 






6 




71 


A = 




~48 




71 


A = 










23vr 


A = 




48 




7T 


A = 






24 



and A is an universal constant which depends only on the type of boundary conditions []TT 

periodic boundary conditions, 
antiperiodic boundary conditions, 
free boundary conditions, 
fixed ++ boundary conditions, 

fixed H — boundary conditions, 

mixed boundary conditions. (27) 

For fixed ++ (or H — ) boundary conditions the spins are fixed to the same (or opposite) val- 
ues on two sides of the strip. The mixed boundary conditions corresponds to free boundary 
conditions on one side of the strip, and fixed boundary conditions on the other. Therefore, 
BK and the mixed boundary conditions are the same in one side of the long strip (fixed to + 

for all spins) and they are different in another side of the long strip (fixed to H 1 — ... for 

BK boundary conditions and free boundary conditions for the mixed boundary conditions). 

Using Kronecker's functions asymptotic form when £ — > and £ — > oo we can obtain 
from Eq. (|TT| ) the free energy per unit length of an infinitely long strip of finite width. In 
the limit £ — > oo (i.e. M. — > oo) for fixed 2M from Eq. ( [P7D one obtains the free energy 
expansion for infinitely long cylinder of circumference 2J\f 

jteo M - 2Uhulk -^ + %\Uf) (2 P )\(2 P + 2) (28) 

V vr 2 J 12 \2J\T) 1440 \2AfJ 
31tt 5 / 1 \ 5 10033tt 7 / 1 \ 7 
~ 24192 V2A/7 ~ 9676800 \2Af) 



This result coincides with that obtained in [14| with the leading finite-size correction to free 
energy — ^(2A/')~ 1 . In the limit £ — > (i.e. A/" — > oo) for fixed we obtain the expansion 
of free energy of infinitely long strip with BK boundary condition of the width Ai 



F TX 

lim —. = Mf bu i k + fi + 



/FZo2M JoaiK JL 24(7W + 1) 



+ E 

P =i 



71 



2(M 
2G 



2p+l . R 



(2p)\(2p + 2) 



(29) 



/ 2C 1 \ 1 
M f--^--ln2) + -ln(l + v^) 

7T 1 

24M* 




2G 


71 1 


71 


f 24 M 


7T 3 \ 




960 J 
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Here the leading finite-width correction to free energy is t^-M 1 . ^From Eqs. (|26"D - (|29|) 
one can see, that 

1. L = 2Af, f = f bulk , f* = 0, A = ~, (30) 

2. L = M, f = f bulk , f* = A, A = ^. (31) 

Our results are consistent with the conformal field theory prediction for the mixed boundary 
condition (see Eq. (§7)) although the mixed boundary condition and the BK boundary 
condition are different in one side of the long strip. 



IV. ASYMPTOTIC EXPANSION OF THE INTERNAL 
ENERGY AND THE SPECIFIC HEAT 

The internal energy per spin and the specific heat per spin can be obtained from the 
partition function Z M ^Af 



U = -2MNTj lnZM ^ = -^2^il lnZM ^ (32) 

1 d 2 
C= 2MAfdT* lnZM ^ 

= MAT \—^diT^ ZM ^-T^ ZM ' m ) ■ (33) 

Let us first consider the internal energy. At the critical point T = T c (fi = 0) the internal 
energy is given by 

U = + In Z 1/2 ,o(0)- (34) 

a [a 

One can note that ^1/2,0 (a 4 ) i s an even function with respect to its argument fx, which imply 
immediately that (^dZi/ 2 ,o(A*)/^A*) Q — 0. Thus we find that internal energy for the finite 

system is equal to its bulk values without any finite-size corrections, namely U = —y2. 
At the critical point T = T c (fi = 0) the specific heat is given by 

AJ\f V2 , r „ /- n1 1 d 2 



Wh,1 + ^ + 2Mv llV *' (35) 



C = -2 - — - 4- tanh 
M M 

The analysis of the Z",^ (0) is a little more involved. Taking the second derivative of Eq. (|TT 
with respect to mass variable 11 and then considering the limit 11 — > 0, we obtain 



^i/2,o(0) = v ! 



= M E ."(^)) cth [M Wo («)] (36) 

A/2,o(UJ „=o 



AT-1 N-l 00 



n=0 n=0 m=l 
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where M = 2{M. + 1), iV = 2j\f, and u}q(x) is the second derivative of ui^{x) with respect to 
fi at criticality 

9 

hi \ 

w oW = - — /7=r^ • 

sin x v 1 + sin x 

Using Taylor's theorem, the asymptotic expansion of the uj'q(x) can be written in the follow- 
ing form 



K2 P X 2 P 



^(2P)I" 

where k 2 = —2/3, k 4 = 172/15, etc. The first sum in Eq. fl3D|) we may transform as 

iV-1 iV-1 ^ c N-l i 

M £ ujp^a) =ME + V E r^TTr (37) 

n=0 n=0 71 n=0 I 

where we have introduce the function f(x) = (jJq{x) — 2/x — 2/(tc — x). This function and all 
its derivatives are integrable over the interval (0, tt). Thus, for the first term in Eq. ( |3"TD we 
may use again the Euler-Maclaurin summation formula, and after a little algebra we obtain 



7V-1 



71=0 71 - 70 p=l * D 



2p 



p(2p)\ 



90 oo T3 1 / 2 1 

where ^f{x) dx = 2 In 2 — 41n7r. The second sum in Eq. ( |3~TD can be written in terms of 
the digamma function ip(x). 

N-l i 

E -—775 = M N + V 2 ) - ^(1/2). ^ 39 ) 

n=0 n "r l l A 

The asymptotic expansion of the digamma function ip(x) is given by (see Appendix [D]) 

oo ol/2 1 

V»(JV + 1/2) = In AT - £(-l)*-£-_. (40) 

p=i ' 

Using the property of the Bernoulli polynomials B^/ 2 , namely, B^+i = 0, Eq. ( [39D can be 
rewritten as 

N-l -I oo T3 1 / 2 1 

V = lniV-y 2p /, -Ml/2). (41) 



Plugging Eqs. ( p8[) and (|4l|) back in Eq. (|37|), we have finally obtain 

ME^) =^{lnAr+iln2-ln7r-V(l/2)} 

n=0 7T L Z J 



^E ^ («) 



p=l 



S J p(2p)\ 
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Let us now consider the second sum in Eq. (p6[). Note that function lo'q(x) can be 
represented as 



1 expiy-^x 2 4 
x P \£l(2p)! J 



(43) 



where coefficients an d «2p are related to each other through relation between moments 
and cumulants (Appendix |B|). Following along the same lines as in the section (3), the 
second sum in Eq. (R3) can be written as 



JV-l oo 



n=0 m=l 

4S {fl 1 AoK)+«l/2)} + (^| + A 2f ^)ln^'« 



7T 



2p 



" 2 ^g P (2p)!V5 



Re Kif°(zAO + 27r£]T 



^(0 

r.l/2 
K 2p -D 2p 



V£Y -1 



! P(2p)! V 5 



where 



Hi/2,0 (0 = -2 In 9^ (0 + Ce + 2 In 2 



(44) 



and Ce is the Euler constant. The differential operators Vl 2p that have appeared here can 
be expressed via coefficients co> 2p = s 2p + ^2 P jx as 

1^4 = &4 -\- 3co>2, 



Substituting Eqs. 



and ([44]) into Eq. (p6|) , we obtain exact asymptotic expansion of 



^i/ 2 ,o(0) 4£ 

^1/2,0(0) VT 



In iV + C s + In 



2 5/2 



7T 



2 In 4 (O 



d 



d 2 



(45) 



p=2 



Plugging Eq. (f45|) back in Eq. fl35|) we finally obtain exact asymptotic expansion of the 
specific heat 



(46) 



where 



11 



8 / , 2 5 / 2 tt\ 4 16. . _ 
Co = - [C E + ln - lnfl 4 £ , 

c 2 = \cx~l + (o 4 2 + el) (i + 4^|) 

2 vr^l^l / 8 3 4 4 21 8 

L/4 = H f7r, U A Or, U A 

£ 270 V 2 2 4 2 4 4 



4* 3 i 2 (mm 7„A /43 « .ft? 



C5 — -C4, 



The ^ expansion of the specific heat has a form 



where 



C7 = -fl + i-V n - M + E^ L ' ( 47 ) 



c = C hi-, 

7T 2 

ci = |c 1 + -(l-ln|), 

2 7T 2 

£ 2 £ 4 
C9 — —Co — — \~ — , 
4 2 tt' 

£ 2 £ 4 

£ 4 £ 2 

Q = 16^ " 2 Cl + to' 

c - 3e V + e V + e r 2 

C5 -"16 C4+ 2" C2 + 2 Cl "5^' 



Typical value of the constants Cq - C3 are given in Table 1, in which the coefficients are 



consistent with those obtained in [05|. 



Using Kronecker's functions asymptotic form when £ — > and £ — > 00 we can obtain 
from Eqs. (^) and Q4"5| ) the specific heat per unit length of an infinitely long strip of finite 
width. In the limit £ —>■ 00 (i.e. .M — > 00) for fixed 2jV the specific heat expansion for 
infinitely long cylinder of circumference 2M can be written as 
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7T ^ V 7T 4/ ^ p(2p)! (2Af) 2p 



7T 7T 



301tt 3 / 1 \ 4 29419vr 5 / 1 \ 6 2759329tt 7 / 1 \° 
10800 V 2Af) ~ 1905120 \2Af) ~ 145152000 \2A? ) ' 

In the limit £ — > (i.e. A/" — ► oo) for fixed Ai we obtain the expression of specific heat of 
infinitely long strip with BK boundary condition of width Ai 

Ti M \ TV J M 

™2k 2p B 2p fir^P- 1 1 

-§(^)^§(^-5) 

8 / 2 3 / 2 
+ - [C E + In + 1 

7T V 7T 



; p(2p)! b; m{m + lyp- 1 



- 



/ 4 7T\/1\ 3 / 2 7T 



V3tt + 18/ \MJ + \3tt + 18 + 21600 














)' 


{m 





(2 7T 43vr 3 \ / 1 
~ 1 5^ + 18 + 7200 J \M 



5 

+ 



Note that the specific heat expansion for infinitely long cylinder contains only even powers 
of M~ l (except of course the leading logarithmic term), while in the specific heat expansion 
for infinitely long strip with BK boundary condition any integer powers of Ai _1 can occur. 

In Fig. 2 we plot the aspect-ratio (£) dependence of the finite-size specific heat correction 
terms Cq, Cx, C2, and C3 for the Ising model with BK boundary condition and those of the 
torus We use the logarithmic scales for the horizontal axis. For large enough £ (3> 1), 



the finite size properties of the Ising model with BK boundary condition and those of the 
torus become the same because the boundaries along the shorter direction determine the 
finite size properties of the system; for both BK boundary condition and the torus, the 
boundary condition along the y axis is the periodic one. 

V. SUMMARY AND DISCUSSION 



In this paper, we have used the method of |20| to derive exact finite-size corrections for 
the free energy F and the specific heat C of the critical ferromagnetic Ising model on the 
Ai x 2J\f square lattice with Brascamp-Kunz (BK) boundary conditions H. We find that 
the finite-size corrections to the free energy and the specific heat are always integer powers 
of A/" -1 (Ai^ 1 ) except of course the leading logarithmic term in the specific heat. In the 
finite-size expansion of the free energy given by Eq. (|18D , only even power of A/" -1 occur, 
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except for the term M . In the finite-size expansion of the specific heat given by Eqs. ( f46|) 
and fl47D , any integer powers of A/" -1 (M.~ r ) can occur. 

We have compared our results with those under toroidal boundary conditions. When 
the ratio £/2 = (M. + l)/2Af is smaller than 1 the behaviors of finite-size corrections for C 
are quite different for BK and toroidal boundary conditions; when ln(£/2) is larger than 3, 
finite-size corrections for C in two boundary conditions approach the same values. In the 
limit M — ► oo we obtain the expansion of the free energy for infinitely long strip with BK 
boundary conditions. Our results are consistent with the conformal field theory prediction 
for the mixed boundary conditions by Cardy |TT] although the definitions of boundary 



conditions in two cases are different in one side of the long strip. It is of interest to know 
under what conditions different boundary conditions could still give the same finite-size 
corrections. 



The results of this paper shows that the method of Ref. pC§ is quite useful for calculating 
exact finite-size corrections for critical systems. It is of interest to apply this method to 
calculate exact finite-size corrections for the Ising model and other free models |20] on various 



lattices with various boundary conditions so that some general features of such finite-size 
corrections could be found. 

This work was supported in part by the National Science Council of the Republic of 
China (Taiwan) under Grant No. NSC 90-21 12-M-001-074. 
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APPENDIX 



APPENDIX A: KRONECKER'S DOUBLE SERIES 



Kronecker's double series can be defined as 122 



K 1/2 '°(r) = V - V 



(m,n)^(0,0) 

In this form, however, they cannot be directly applied to our analysis. We need to cast 
them in a different form. To this end, let us separate from the double series a subseries with 
m = 

K iA0 M = P l T ^1 _ P l T T e-™ 

{-2m)P^ Q nP (-27ri)P^ ^ z (n + rm)p- 

Here the first sum gives nothing but Fourier representation of Bernoulli polynomials 

B" = - P \ Y- . (Al) 

The second sum can be rearranged with the help of the identity 



p\ 



Y r = P Y(n + l /2Y^e 2mz ^ +1 ^ , 



h2m)P^ z (z + n)P n=Q 
which can easily be derived from following equation 

p 2mza co -r +00 — 2irina 
JL = _y e 2^(n+a) = _L y e _ (A2 ) 

e^-1 ^ 2« n f^ 00 z + n 

by differentiating it p times. The final result of our resummation of double Kronecker sum 
is 

00 

KJ/ 2 '°(t) = B l J 2 - P YY(n + 1/2)^ e 2mmr{n+1/2) . 

Considering the Kronecker sums with pure imaginary aspect ratio, r = z£, we can further 
rearrange this expression to get summation only over positive m > 1 

oo oo 

K 2 - K 2p 2 '°(^) = 4p Y Z> + 1/2)*" 1 e~ 2 ^ n+1 l 2 \ (A3) 

m=l n=0 



APPENDIX B: RELATION BETWEEN MOMENTS AND CUMULANTS 

Moments Z2k and cumulants which enters the expansion of exponent 

f 00 ^2k ~\ oo g,2k 
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are related to each other as 
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%2 — F 2 , 

Z4 = F4 + 3-F 2 , 

Z 6 = F 6 + 15F 2 F 4 + 15F|, 

Z 8 = F 8 + 28F 2 F 6 + 35F 4 2 + 210F 2 2 F 4 + 105F 2 4 , 

where summation is over all positive numbers and different positive numbers 

{k\, . . . , k r } such that fciii + . . . + k r i r = k. 



APPENDIX C: REDUCTION OF KRONECKER'S DOUBLE SERIES TO THETA 

FUNCTIONS 



Let us consider Laurent expansion of the Weierstrass function 



(n,m)^(0,0) 



1 



-2 + JL a p( 



T)Z 



(z — n — rm) 2 (n + rm) 2 

2p-2 



p=2 



The coefficients a p (r) of the expansion can all be written in terms of the elliptic ^-functions 
with the help of the recursion relation pB 



p (p-3)(2p + l) 
where first terms of the sequence are 



(a 2 a p _2 + a 3 a p _ 3 + . . . + a p _ 2 a 2 ) 



a 2 = i(^f-^! + M 



3^47' 



■I J; 



(I I 



3 U 2) 



05 — Yil^CtsJ, 

a 6 = ±(2al + 3a 2 3 ), 



Kronecker functions K 2 p°(r) are related directly to the coefficients a p (r) 



K 2 P (t) 



(2p)! a P (r) 

;-4vr 2 )p (2p - r 



Kronecker functions K 2 p 2 '°(r) can in their turn be related to the function K 2 p°(r) by means 
of simple resummation of Kronecker's double series 
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K|'°(r)=2 1 ^Kj (r/2)-K^(r). 

Thus, Kronecker functions K^ 2 ' ^) can all be expressed in terms of the elliptic ^-functions 
only. For practical calculations the following identities are also helpful 



26 2 2 {2t) 




6l(r/2) 


= 28 2 s , 


2#a (2r) 


= 01 + 91 


6l(r/2) 


= 6 2 + 3 , 


20|(2r) 


= 29 3 9i, 


0l(r/2) 


= 61- 61 



^From the formulas above we can easily write down the Kronecker functions that have 
appeared in our asymptotic expansions 

K I'°(r) = ^(|^-«, 



4'°(t) = -±(e$ + el)(%el + ei9i). 

Note that when £ — > 00 we have limits 62 — * 0, 64 — > 1, 63 — > 1 . The case £ — > can 
be obtained by using Jacobi's imaginary transformation of the 6 - functions. In this case 
60 — > -k=, 6a — > and 63 — > -K= and the Kronecker's function can again be reduce to the 

Bernoulli polynomials. 



APPENDIX D: ASYMPTOTIC EXPANSION OF THE DIGAMMA FUNCTION 

ip(N + a) 

Let us start with well known expansion of the digamma function i/j(N) 

ip(x) =hix - ^- ~E >2 '' 



2x fz^ 2p x 2p 

00 R 1 

= lnx-E(-l) P ^4- ™ 

p= i V % p 

Plugging in the above expansion x = N + a and expand the resulting factors In (1 + a/N), 
(1 + a/N)~ p in powers of iV _1 we obtain 

00 a p 00 00 , (n + k — lV o k 

00 o p 00 - (I — IV ry l ~ p 

= iniv - E(-i) p ^ - E E(-i)'Bp7r-in V 

p^i p Np £ip^i (i -py.pi N l 

= \nN-Y Y(-l)% } 7', . (D2) 
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Using the relation between Bernoulli polynomials B° and Bernoulli numbers B p 



we finally obtain Eq. 



Br= S B '(r4p a, " P ' (D3) 



OO gQ -| 

^(N + a)=\nN- £(-1)* ^ Wp - (D4) 
p=i P iV 
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FIGURE CAPTIONS 



FIG. 1. The Brascamp - Kunz boundary conditions for M x 2N lattice. Here M = 7 
and 2N = 8. 

FIG. 2. Aspect-ratio (£) dependence of finite-size correction terms for the specific heat 
of the square lattice Ising model with Brascamp-Kunz boundary conditions (solid lines) and 
toroidal boundary conditions (dashed lines): (a) Co, (b) Ci, (c) C 2 , and (d) C 3 . 



21 



2N 

Fig 1 . Braskamp-Kunz boundary conditions. M x 2N = 7 x 8 



Table 1. Values of the coefficients c - c 3 for various values of the ra- 
tio i={M + 1)1 M. Here are presented (c - 2)/8 and a/8 (i = 1,2,3) 
for the convenience of comparison with results of Janke and Kenna [cond- 
mat/0103332 (2001)], in which 1/ p is corresponding to £ in the present paper. 





1/2 


1 


2 


(co - 2)/8 


-0.349 694 206 9... 


-0.350 879 733 2... 


-0.376 674 233 4... 


ci/8 


0.291 838 983 9... 


0.290 653 457 6... 


0.264 858 957 4... 


c 2 /8 


0.180 950 438 7... 


0.175 784 345 6... 


0.125 896 137 8... 


c 3 /8 


-0.074 847 143 3... 


-0.069 681 050 2... 


-0.019 792 842 4... 
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